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Abstract. We derive tropospheric column BrO during the ARCTAS and ARCPAC field cam- 
paigns in spring 2008 using retrievals of total column BrO from the satellite UY nadir sensors 
OMI and GOME-2 using a radiative transfer model and stratospheric column BrO from a 
photochemical simulation. We conduct a comprehensive comparison of satellite-derived tro- 
5 pospheric BrO column to aircraft in-si tu observations of BrO and related species. The aircraft 
profiles reveal that tropospheric BrO, when present during April 2008, was distributed over 
a broad range of altitudes rather than being confined to the planetary boundary layer (PBL). 
Perturbations to the total column resulting from tropospheric BrO are the same magnitude 
as perturbations due to longitudinal variations in the stratospheric component, so proper ac- 
1 0 counting of the stratospheric signal is essential for accurate determination of satellite-derived 
tropospheric BrO. We find reasonably good agreement between satellite-derived tropospheric 
BrO and columns found using aircraft in-situ BrO profiles, particularly when satellite radi- 
ances were obtained over bright surfaces (albedo > 0.7), for solar zenith angle < 80° and clear 
sky conditions. The rapid activation of BrO due to surface processes (the bromine explosion) is 
1 5 apparent in both the OMI and GOME-2 based tropospheric columns. The wide orbital swath 
of OMI allows examination of the evolution of tropospheric BrO on about hourly time inter- 
vals near the pole. Low surface pressure, strong wind, and high PBL height are associated 
with an observed BrO activation event, supporting the notion of bromine activation by high 
winds over snow. 

20 1 Introduction 

Bromine plays an important role in tropospheric ozone chemistry and the resulting oxidation ca- 
pacity of the polar boundary layer. Bromine radicals catalytically destroy ozone, leading to nearly 
complete removal near the surface that is termed an ozone (O 3 ) depletion event (ODE). Once O 3 
is depleted, high levels of reactive halogen species including atomic bromine (e.g. Br) become the 
25 primary oxidants for many species, including methane (CH 4 ) and mercury (Hg) (e.g. Simpson et aL, 
2007b; Schroeder et aL, 1998). During mercury depletion events (MDEs), Br and BrO are thought 
to oxidize elemental mercury to more reactive gaseous mercury that deposit to the polar ecosystem 
(Schroeder et aL, 1998; Lu et aL, 2001 ; Ariya et aL, 2004; Douglas et aL, 2005; Holmes et aL, 2010). 

A primary source of active bromine in the polar boundary layer is thought to be bromide (Br~) 
30 in the condensed phase, which is transformed to gaseous molecular bromine (Br 2 ) by heterogenous 
chemistry. Potential sources of bromide are sea salt aerosols (Fan and Jacob, 1992), surfaces of first 
year sea ice (Simpson et aL, 2007a; Wagner et aL, 2007) or newly formed sea ice (Jones et aL, 2006), 
frost flowers (Kaleschke et aL, 2004), or blowing snow triggered by strong winds (Yang et aL, 2008, 
2010; Jones et aL, 2009; A. E. Jones et aL, 2010; Begoin et aL, 2010; Theys et aL, 2011; Toyota et 
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35 al., 201 1 ). Br 2 released to the atmosphere is rapidly photolyzed to yield bromine radicals (Br), i.e. 


HOBr + Br“ +H + ?=± H 2 OT*Br 2 (condensed phase) ( 1 ) 

Br 2 + h i/ 2 Br (gas phase). ( 2 ) 

The bromine radicals then lead to catalytic ozone loss, i.e. 

Br + 0 3 ~>Br0~h0 2 (3) 

40 BrO + XO -» BrX + 0 2 (X = Br, Cl, I, and OH) (4) 

BrX + h^— >Br + X (5) 


Barrie et al. (1988) found that filterable bromine is strongly anti-correlated with the abundance 
of O 3 in the polar boundary layer. Following this work, many additional measurements of bromine 
species along with O 3 have been made in polar environments (e.g. Bottenheim et al., 1990; Barrie 
45 et al., 1994; Hausmann and Platt, 1994; Tuckermann et al., 1997; Martinez et al., 1999; Ridley et 
al., 2003). The association of elevated BrO and depleted surface O 3 during Arctic spring is well 
established (e.g. Platt and Honninger, 2003, and references therein). 

Bromine monoxide (BrO), an intermediate in the catalytic loss of ozone, is the most commonly 
observed active bromine species. BrO absorbs ultraviolet (UV) radiation, which enables measure- 
50 ment using remote sensing techniques. Reported in-situ and remotely-sensed ground-, balloon-, 
and aircraft-based BrO measurements include: ( 1 ) boundary layer mixing ratio with Long Path- 
Differential Optical Absorption Spectroscopy (LP-DOAS) (Hausmann and Platt, 1994; Tuckermann 
et al., 1997; Martinez et al., 1999); (2) mixing ratio by aircraft in-situ Chemical-Ionization Mass 
Spectrometry (CIMS) (Neuman et al., 2010; Liao et al., 201 la,b); (3) boundary layer and free tropo- 
55 spheric BrO column using DOAS on aircraft (Me Elroy et al., 1999; Prados-Roman et al., 2011); (4) 
total, stratospheric, and tropospheric columns using ground-based zenith sky and direct sun DOAS 
(Friefi et al., 1999; Schofield et al., 2004, 2006; Hendrick et al., 2007, 2008; Theys et al., 2007); (5) 
boundary-layer column with ground-based MAX-DO AS (Honninger, 2004; Simpson et al., 2007a; 
Donohoue et al., 2010; Friefi et al., 2011); ( 6 ) balloon profiles using a variety of spectroscopic 
60 methods (Fitzenberger et al, 2000; Pfeilsticker et al., 2000; Pundt et al., 2002; Dorf et al., 2008). 
Despite the many measurements of BrO and related species obtained by various techniques, sig- 
nificant uncertainties remain regarding the importance of very short lived source compounds on the 
stratospheric bromine budget (e.g. Sect. 1.3.3. 3 of WMO, 2011) as well as the magnitude of the 
global, ubiquitous, background level of tropospheric BrO (e.g. Sect. 5.3 of Theys et al., 2011). 

65 Space-based observation of BrO offers an excellent tool for studying Arctic polar bromine chem- 
istry. Satellite observations provide global coverage, far superior to the spatial coverage available 
from ground- and aircraft- based measurements. Limb-sounding observations from the Microwave 
Limb Sounder (MLS) (Livesey et al, 2006; Kovalenko et al., 2007) provide quantification of BrO 
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profiles in the upper and middle stratosphere. The MLS observations of BrO imply a significant 
70 contribution to stratospheric bromine from sources other than long-lived CH 3 Br and halons. Limb 
observations from the Scanning Image Absorption SpectroMeter for Atmospheric CHartographY 
(SCI AM ACHY) (Sinnhuber et al., 2005; Sioris et al., 2006; Rozanov et al., 2011) extend into the 
lowermost stratosphere (LMS) and nadir measurements from SCIAMACHY constrain total column 
BrO, including contributions from the stratosphere and troposphere. The retrievals of SCIAMACHY 
75 BrO by Rozanov et al. (201 1) and Sinnhuber et al. (2005) imply a limited role for supply of bromine 
by very short lived (VSL) compounds to the LMS and a considerable burden of global, ubiquitous, 
background BrO whereas the retrieval of SCIAMACHY BrO described by Sioris et al. (2006) im- 
plies a larger role for supply of stratospheric inorganic bromine, Br y , by VSL compounds and a 
much smaller level for background tropospheric BrO (e.g. Sect. 2.5.2. 1 of WMO, 2007). 

80 Nadir-viewing instruments on polar-orbiting satellites provide multiple daily observations of total 
column BrO at high latitude. The high surface albedo of polar regions provides good sensitivity 
to tropospheric BrO, including that near the surface (Wagner and Platt, 1998; Theys et al., 201 1). 
Chance (1998) and Hegels et al. (1998) retrieved BrO total vertical column densities (VCD) from 
nadir radiances in the UV obtained by Global Ozone Monitoring Experiment (GOME) (GOME 
85 Users Manual, 1995) on the European Space Agency (ESA) European Remote Sensing 2 (ERS- 
2) satellite. The early GOME observations showed large enhancements over Hudson Bay during 
spring 1997, which was attributed to bromine release from the surface (Chance, 1998). Estimates of 
BrO total column amount have been subsequently derived from other nadir-viewing satellite sensors 
including the Ozone Monitoring Instrument (OMI) (Levelt et al., 2006) on the US National Aeronau- 
90 tics and Space Administration (NASA) Aura satellite, SCIAMACHY (Bovensmann et al., 1999) on 
the ESA Environmental Satellite (EnviSat), and the second Global Ozone Monitoring Experiment 
(GOME-2) instruments (Munro et al., 2006) that are flying on the series of European Meteoro- 
logical Operational Satellites (EuMetSat MetOp). Retrievals of VCD BrO from different satellite 
instruments and by various groups tend to be in fairly close agreement. All satellite instruments 
95 report levels considerably larger than predicted by models that account only for stratospheric supply 
from the decomposition of methyl bromide and halons. The focus of the scientific community since 
publication of early GOME observations has been quantifying the degree to which this “excess BrO” 
resides in the troposphere, the stratosphere, or perhaps both regions of the atmosphere. 

Tropospheric BrO columns can be derived from satellite observations using the residual method, 
1 00 in which an estimate of the stratospheric BrO column is subtracted from the satellite-derived VCD of 
BrO. Wagner and Platt (1998) and Richter et al. (1998) estimated tropospheric VCD of BrO using 
this approach from GOME observations, based on simple assumptions for stratospheric BrO. These 
studies and many others assumed that stratospheric BrO is zonally symmetric and that bromine was 
supplied to the stratosphere solely by the decomposition of long-lived organic compounds. Theys et 
105 al. (2009) and Salawitch et al. (2010) showed that the stratospheric VCD of BrO can exhibit strong 
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gradients with respect to longitude at high latitude during spring and Salawitch et al. (2010) ques- 
tioned prior estimates of residual tropospheric BrO found assuming that the stratospheric burden 
was zonally symmetric. Theys et al. (201 1) derived tropospheric BrO columns from GOME-2 spec- 
tra with a model-based stratospheric BrO climatology (Theys et al., 2009) similar to that used here 
1 1 0 and described below. 

Many studies related to bromine and ozone chemistry have been conducted using satellite-derived 
tropospheric BrO columns. Wagner and Platt (1998) reported elevated regions of BrO vertical col- 
umn density in the Arctic and Antarctic regions. They noted these enhancements were likely due to 
increased abundance of tropospheric BrO, rather than a stratospheric disturbance, based on a vari- 
115 ety of factors including the correlation between enhanced columns of BrO and the O 2 -O 2 collision 
complex. Wagner et al. (2001) showed that elevated BrO column amounts observed by GOME 
were correlated with low ozone in the boundary layer observed in-situ at Ny-Alesund (Spitsber- 
gen), Norway. In other studies, spatial and temporal features of ODEs have been simulated using 
3 -dimensional regional chemical transport models and GOME-derived tropospheric BrO columns 
1 20 (Zeng et al., 2003, 2006). Connections between Br O-rich air masses and first-year sea ice have been 
indicated with back-trajectory analyses using SCIAMACHY data (Wagner et al., 2007). Transport 
of a large BrO plume near the North Pole is also reported by Begoin et al. (2010). A back trajectory 
study using satellite-derived tropospheric BrO columns indicated that ODEs can be differentiated 
into locally activated and transport driven events (Koo et al, 2011). 

1 25 Despite the numerous studies of tropospheric polar bromine chemistry using satellite BrO obser- 
vations, estimation and interpretation of tropospheric BrO information from space presents ongoing 
challenges. To properly estimate tropospheric BrO column amounts, the stratospheric contribution 
to the satellite-derived total column must be accurately represented (e.g. Theys et al., 2009; Salaw- 
itch et al., 2010). The global, ubiquitous background tropospheric level of BrO inferred from the 
1 30 satellite record is sensitive to the amount of Br y delivered to the stratosphere by VSL bromocarbons 
(Salawitch et al., 2005). Furthermore, low solar elevation angles in the early polar spring lead to 
large uncertainties in satellite total BrO column retrievals (see below). The presence of clouds fur- 
ther complicates the retrieval of tropospheric BrO from satellite observations. Theys et al. (201 1) 
have recently addressed many of these issues. They showed maps of tropospheric BrO columns 
135 derived using a method similar to that described below, and evaluated these columns using ground- 
based measurements of the tropospheric and stratospheric contributions to the total column. Our 
study builds upon the work of Theys et al. (201 1) by further exploring the issues that affect quan- 
tification of tropospheric BrO columns and by using aircraft measurements of BrO to evaluate the 
satellite-derived tropospheric columns. 

140 Salawitch et al. (2010) provided an initial analysis of BrO observations obtained by instruments 
aboard the NASA and NOAA aircraft (Neuman et al., 2010; Liao et al., 201 lb) and the OMI satel- 
lite instrument (Kurosu and Chance, 2011) during spring 2008. They simulated the stratospheric 
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BrO column assuming VSL bromocarbons supply between 5 and lOppt of Br y to the stratosphere 
in the form of product gas injection (PGI), in addition to the ppt of Br y supplied by source gas 
145 injection (SGI) of the VSL species CH 2 Br 2 . On the other hand, Theys et al. (2011) simulated the 
stratospheric burden of BrO assuming 1 ppt of PGI and 5 ppt of SGI (all represented as CH 2 Br 2 ). 
The stratospheric BrO burden found by Salawitch et al. (2010) for their “best case” simulation of 
7 ppt from PGI and 2 ppt from SGI (see their Figure S7) is considerably larger (^27 % overall dif- 
ference, with quite a bit of geographic variability) than the stratospheric BrO burden used by Theys 
1 50 et al. (201 1) as shown in Figure SI. Although we do not compare to independent observations 
of BrO here, such comparisons to aircraft and balloon-borne observations of BrO in Salaw- 
itch et al. (2005) support the use of our values for the delivery of stratospheric bromine by PGI 
and SGI. Other observations of BrO, such as some of those discussed in Section 2.5.2 of WMO 
(2007), support smaller values for these parameters. 

1 55 Salawitch et al. (2010) concluded that high column amounts of BrO derived from satellite obser- 
vations could, in some cases, be attributed to compression of stratospheric air to low altitudes (high 
pressures) and that this condition, not previously recognized in the analysis of the satellite BrO 
record, should be considered to properly relate the satellite record to surface events. The magnitude 
of the stratospheric enhancement is sensitive to the supply of stratospheric Br y from VSL bromocar- 
1 60 bons as well as the pathways of this transport (i.e, PGI versus SGI; see Chapter 2 of WMO (2003) 
for a detailed discussion of these pathways). Salawitch et al. (2010) showed agreement, to within 
uncertainties, between total column BrO measured by OMI and the sum of modeled stratospheric 
and aircraft-measured tropospheric BrO partial columns. However, their treatment of the tropo- 
spheric column did not explicitly account for tropospheric air mass factors (AMFs) and thus must 
1 65 be treated with caution. Our study builds on this prior work by using a radiative transfer model to 
calculate tropospheric AMFs and also by considering the effects of clouds, surface reflectivity, and 
viewing geometry on the evaluation of the bromine budget for Arctic spring 2008. 

Validation of satellite -derived BrO columns and their relationship to ODEs was a goal of two 
Arctic field campaigns conducted in 2008 under the auspices of the International Polar Year (IPY). 
1 70 The Arctic Research of the Composition of the Troposphere from Aircraft and Satellite (ARCTAS) 
mission was conducted in April and June-July 2008 by NASA (D. J. Jacob et al., 2010) and the 
US National Oceanic and Atmospheric Administration (NOAA) Aerosol, Radiation, and Cloud Pro- 
cesses affecting Arctic Climate (ARCPAC) mission took place in April 2008 (Brock et al., 2011). 
Various in-situ aircraft measurements of trace chemicals, including Os and reactive bromine species 
175 (BrO, Br 2 , and soluble bromide), were made during these campaigns. Neuman et al. (2010) and 
Liao et al. (2011a) provide an overview of the in-situ measurements of BrO and related species 
obtained during ARCTAS and ARCPAC. 
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2 Data description 


We use a variety of aircraft and satellite measurements and model outputs to calculate and evaluate 
1 BO satellite-derived tropospheric BrO columns. Tropospheric BrO columns are inferred from aircraft 
measurements of in -situ BrO (henceforth referred to as “in-situ columns”, as described in Sect. 2.1). 
We derive tropospheric vertical column densities (VCDs) using BrO slant column densities (SCD) 
retrieved from OMI and GOME-2 (see Sect. 2.2) and a model simulation of the stratospheric BrO 
column as detailed in Sect. 2.3. OMI rotational Raman scene pressures and MODIS cloud products 
1 85 (see Sect. 2.4) are used to assess cloud effects on the derived tropospheric BrO columns. Data from 
a global reanalysis, described in Sect. 2.5, are used to examine relationships between observed BrO 
enhancement events and meteorological conditions. 

2,1 Aircraft in-situ measurements 

Aircraft in-situ measurements of BrO and various other trace gas concentrations were made from 
190 the NASA DC-8 aircraft during ARCTAS and the NOAA WP-3D aircraft during ARCPAC. An 
overview of the instruments, flights, and mission goals and accomplishments are provided in the 
ARCTAS (D. J. Jacob et al., 2010) and ARCPAC (Brock et al., 201 1) overview papers. BrO was 
measured using Chemical Ionization Mass Spectrometer (CIMS) instruments during both experi- 
ments (Neuman et al., 2010; Liao et al., 201 1 b). 

1 95 DC-8 CIMS measurements of BrO were made once every 30 s, and WP-3D CIMS measurements 
were acquired once every 2 s. Measurement uncertainties for BrO are ±40 % with a detection limit 
of 3 pptv for WP-3D data and ±40 % and a detection limit of 2-5 pptv for DC-8 data (Neuman et 
al., 2010; Liao et al., 201 lb). In this study, we use 1 min averaged data for both DC-8 and WP-3D 
measurements. The spatial resolution of the 1 min averaged aircraft data is approximately 10 km. 
200 We estimate tropospheric BrO columns from the aircraft BrO mixing ratio measurements for each 
ascent and descent. Figure 1 shows the flight tracks and locations of the 1 6 (29) profiles made from 
the DC-8 (WP-3D). BrO profiles are available for 16 and 17 April 2008 from the DC-8 and 12, 15, 
18, 19 and 21 April 2008 from the WP-3D. We also examine measurements of other quantities 
for DC-8 flights of 5 and 8 April 2008, for which BrO was below the instrument detection limit. 
205 To calculate tropospheric columns from the aircraft data, the profiles were binned on an altitude grid 
with 500 m resolution from the surface to 7.5 km (BrO measurements are available only up to 7.5 km 
altitude). For each bin, the median value was selected as the BrO mixing ratio at that altitude. If the 
median value was less than 0, then we assumed the mixing ratio was 0. 

This sampling strategy was chosen to properly represent the signature of the ARCTAS and 
210 ARCPAC BrO profiles, which show that the tropospheric column is domimated by contribu- 
tions from above the top of the convective boundary layer (i.e. 9 the region of constant potential 
temperature; see Neuman et al. (2010) and Salawitch et a!. (2010)) and, when BrO enhance- 
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ments occur, they are present in relatively thick altitude segments (i.e., extent larger than 500 
m). Our results are insensitive to use of mean BrO within the layers rather than median BrO. 

21 5 Figure 2a shows the entire collection of DC-8 BrO measurements and the median DC-8 (hence- 
forth referred to as the composite) profile. Figure 2b and 2c shows samples of different types of 
DC-8 BrO profiles collected during April 2008. Among the 16 DC-8 BrO profiles, 8 have shapes 
similar to the DC-8 composite profile, while the others show elevated BrO near the surface and/or 
in layers at altitudes from 2 to 4 km. WP-3D measurements of BrO exhibit larger variability than 
220 those from the DC-8 (Fig. 2d, e, and f). 

The DC-8 and WP-3D instruments reported much lower mixing ratios of BrO near the sur- 
face than have been measured in the past by LP-DOAS instruments in the springtime Arctic 
boundary layer (Hausmann and Platt, 1994; lYickermann et aL, 1997; Martinez et al., 1999). 
The median value of surface BrO at Barrow, Alaska (71° N, 156° W) during spring 2009 was 
225 measured to be ^4 ppt by a CIMS instrument and ppt by a DOAS instrument (Figure S2). 
We have no reason to suspect that conditions during spring 2008 were appreciably different 
than conditions during spring 2009. The low values of BrO observed during ARCTAS and 
ARCPAC could be due to preferential sampling of 0 3 depleted air by the two aircraft (which 
would titrate inorganic bromine from BrO to Br) or the fact the aircraft generally sampled 
230 to an altitude of 100 m above the surface (Neuman et al., 2010). While it is tempting to argue 
that perturbations to the BrO profile due to surface release are confined to the lowest region 
of the marine boundary layer (a region inaccessible to the two aircraft), ground-based MAX- 
DOAS observations from Barrow, Alaska obtained during April 2008 suggest that when BrO 
was elevated, the perturbation extended to ^1 km altitude (Salawitch et al., 2010). In contrast, 
235 MAX-DOAS near Barrow obtained during April 2009 by another group suggest the enhance- 
ment to BrO usually occurs in the lowest 300 m, with one case where elevated BrO extended 
to 500 m altitude (Friefi et al., 2011). A retrieval of the BrO profile using radiances mea- 
sured in the Arctic during April 2007 by an airborne limb scanning mini-DOAS instrument 
also suggest elevated BrO is confined to the boundary layer (Prados-Roman et al., 2011). The 
240 ARCTAS and ARCPAC in-situ profiles of BrO show, consistent with MAX-DOAS data from 
April 2008, that the highest mixing ratios of BrO tended to be observed above the boundary 
layer. Nonetheless, the profiles of BrO shown in Fig. 2 almost certainly reflect admixtures 
of air that has been influenced by recent surface release of active bromine with air that has 
achieved a background level of BrO, reflecting dilution of the surface signal. A plausible ex- 
245 planation for the tendency for aircraft BrO to be lower than surface BrO is the influence of 
mixing. Satellite measurements of total column BrO will be influenced by this mixing process, 
as well. Below, we quantify the impact of elevated surface BrO on the total column. 

To estimate the tropospheric BrO column from aircraft measurements, we must first make as- 
sumptions about mixing ratios between the surface and the lowest altitude sampled by the aircraft. 
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250 We only take profiles with aircraft minimum altitudes less than 500 m. Then, we assume that the 
BrO mixing ratio in the lowest bin (surface to 500 m) is the median BrO mixing ratio between the 
lowest aircraft altitude and 500 m. When the aircraft did not sample up to the tropopause, we made 
assumptions about mixing ratios between the highest aircraft altitude and the tropopause. Here, we 
use the upper part of the DC-8 composite profile to fill empty upper bins. However, the DC-8 corn- 
255 posite profile only goes up to 7.5 km, the highest altitude where BrO is sampled by the aircraft. We 
assume that the BrO mixing ratio between 7.5 km and the tropopause is zero. Aircraft data suggest 
that BrO mixing ratios at these altitudes are very small for the sampled air masses (see Fig. 2). It 
is possible that stratospheric to tropospheric transport of air mass could supply BrO to the upper 
troposphere, particularly along the western flank of Arctic low pressure systems (i.e. after low alti- 
260 tude tropopause systems pass over a region) (Salawitch et al., 2010). This idea is speculative and is 
not considered below. Finally, layer column BrO amounts (BrO column amounts for vertical bins) 
are estimated using pressures and temperatures from aircraft measurements and integrated from the 
surface to the tropopause, to provide a tropospheric column BrO abundance. We use the tropopause 
height from the Modem Era Retrospective- Analysis for Research and Application (MERRA) data 
265 set (see Sect. 2.5). 

Tables 1, 2a, and 2b provide information related to the aircraft BrO profiles including the bottom 
and top altitudes of the aircraft in-situ profiles (Alts and Alfr, respectively), the MERRA tropopause 
height (Htrp), and the in-situ tropospheric BrO column (Col. in-situ) obtained by integrating the 
aircraft profiles. Tables 3, 4a, and 4b provide solar zenith angle (SZA), OMI reflectivity at 331 nm 
270 (7^331 nmX the difference between terrain and OMI cloud pressure (A P c , explained in Sect. 2.4), 

the satellite-derived tropospheric columns from OMI (CoI.omi) and GOME-2 (CoI.gome-2) cor- 
responding to the in-situ aircraft BrO profiles, and the ratios of satellite-derived tropospheric BrO 
columns to the in-situ BrO columns (RatiooMi and RatiOGOME-2) (see Section 2.2 and 3.2 for a 
description of the satellite data and related parameters in these tables). 

275 The lowest altitudes sampled by the DC-8 and WP-3D during their descents over BrO- 
enhanced regions were 75 and 77 m, respectively. In our analysis we use the composite DC-8 
or WP-3D profile to extrapolate between the lowest sampled altitude and the surface, for each 
formulation of in-situ column BrO shown in Tables 1, 2a, and 2b. At times, surface BrO can 
reach mixing ratios as high as 40 pptv (Liao et al., 2011a). We have assessed the impact of ele- 
280 vated surface BrO on our analysis of aircraft, satellite, and modeled stratospheric columns by 
conducting a probability distribution function for daytime surface BrO, observed at Barrow, 
Alaska. Two thirds of the time, surface BrO is below 8 pptv (Figure S2). A uniform distribu- 
tion of 8 pptv of BrO between the surface and 75 m altitude would contribute 0.18xl0 13 cm" 2 
to the column, an amount much smaller than the in situ and satellite-based columns discussed 
285 throughout the paper (see caption, Figure S2). Levels of BrO reaching 40 pptv below the air- 
craft would contribute IxlO 13 cm " 2 to the column if the BrO where present, at this amount, 


10 



uniformly between the surface and 75 m altitude. The PDF analysis shows that while surface 
BrO did reach 40 ppt in spring 2009, such occurences were rare. While layers of highly el- 
evated BrO below the aircraft could on occaison compromise our comparisons, our overall 
290 conclusions are robust because surface measurements indicate only on rare occaison are BrO 
enhancments large enough to significantly perturb the column. 

In addition to BrO, other aircraft trace gas measurements are useful for inferring halogen chem- 
istry and air mass characteristics. For example, Bi’2 was measured by CIMS instruments from both 
aircraft. Laboratory studies revealed that the Br2 signals include contributions from IiOBr that re- 
295 acted on inlet surfaces. Therefore, Br 2 represents the lower limit of HOBr + Br2 (Neuman et al„ 
2010). Henceforth, we refer to this measurement as “active bromine”. Soluble bromide was also 
measured from the DC-8 using mist chamber/ion chromatography. Soluble bromide may include 
inorganic bromine compounds such as HBr, HOBr, particulate bromine, and possibly BrO and Br2 
(Ridley et al., 2003; Dibb et al., 2010; Neuman et al., 2010; Liao et aL, 201 lb). Measurements of O3 
300 mixing ratios by chemiluminescence were made from both aircraft. CO was measured using tunable 
diode laser absorption spectroscopy on the DC-8 and UV fluorescence on the WP-3D (D. J. Jacob et 
aL, 2010; Brock et al.,2011). 

2.2 Total BrO slant column density from OMI and GOME-2 

OMI is a nadir-viewing ultraviolet and visible (UV/Vis) sensor (Levelt et al., 2006) aboard the NASA 
305 Aura satellite that is in a sun-synchronous orbit with an overpass of 01:38 p.m. local time. The 
spectral resolution in the OMI UV-2 channel used to retrieve BrO columns is approximately 0.5 nm. 
The OMI swath width is about 2600 km. The pixel size of OMI UV-2 channel is approximately 
1 3 x 24 km 2 at the swath center and significantly larger at the swath edges. With its wide swath, OMI 
provides multiple daily observations at high latitudes in spring and daily global coverage at low and 
31 0 middle latitudes. An obstruction outside the instrument that produces radiance errors (known as the 
“row anomaly”) reduced the swath coverage mainly after May 2008; it does not significantly affect 
the observations shown here (Claas et al., 2010). 

OMI BrO SCDs are retrieved by directly fitting backscattered UV radiances to absorption cross- 
sections of BrO (the target gas), NO2, HCHO, and SO2 as well as inelastic rotational-Raman scat- 
31 5 tering (also known as the Ring effect) using a non-linear least-squares approach (Chance, 1998). 
The spectral fitting window for the OMI algorithm is 319 to 347.5 nm. BrO cross sections from 
Wilmouth et al. (1999) are used. The total fitting uncertainty of OMI BrO total column typically 
ranges from 15 to 51 % (Salawitch et al., 2010; Kurosu and Chance, 201 1). 

GOME-2 is a 4 channel UV/Vis nadir viewing instrument operating on Metop-A since January 
320 2007 (Callies et al., 2000). It has a local equator crossing time of 9:30 am in the descending node, 

a swath width of 1920 km, and a spatial resolution of 40 x 80 km 2 . For the BrO SCD retrieval used 
here (Begoin et al., 2010), measurements in the window 336 to 347 nm are used, where GOME-2 
has a spectral resolution of about 0.3 nm. SCDs of BrO are retrieved using the standard DOAS 
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approach. Absorption due to BrO (Wahner et al., 1988), O 3 , NO 2 , and the effects of rotational- 
325 Raman scattering are included. The uncertainty of the GOME-2 BrO total SCDs is 10 to 30% 
depending on solar zenith angle (SZA) and surface albedo and the uncertainty has both random and 
systematic contributions from spectral interferences and the cross-sections. As a result of throughput 
loss of the UV channels, random errors have increased since launch. However, this effect is not 
significant during the time period studied here (Dikty et al., 201 1 ). 

330 All of the analyses conducted in this study are based on the assumption that the magni- 
tude of satellite-derived slant column is correct. A detailed analysis of the errors involved in 
deriving satellite-derived slant columns of BrO is beyond the scope of this paper. There are 
potentially large uncertainties associated with the derivation of both slant and vertical absolute 
column BrO related to the choice of spectral fitting windows and various DOAS parameters 
335 (e.g., polynomial order, orthogonalizations) used in retrievals as well as assumed BrO vertical 

profiles (G. Mount, private communication, 2011; DLR, 2009). 

For all analyses conducted in this study, we use level 2 (time-ordered) satellite data. Level 2 data 
provide the best temporal match between satellite observations and aircraft in-situ and ground-based 
observations. 

340 2.3 Stratospheric BrO column 

A model simulation of stratospheric bromine species is used to estimate the spatial structure of the 
stratospheric BrO column. Our approach is similar to that described by Salawitch et al. (2010), but 
differs in that they showed results for a fixed local solar time of 01 :30 p.m. whereas here we calculate 
BrO along each OMI or GOME-2 orbit, which spans a range of local solar times near the pole. A 
345 photochemical steady state (PSS) model is constrained to profiles of temperature, O 3 , NO y , Cl y , 
H 2 O, CH 4 , sulfate aerosol surface area, etc. output from a run of Whole Atmosphere Community 
Climate Model (WACCM) (Garcia et al., 2007) conducted using meteorological fields for spring 
2008. Salawitch et al. (2010) provide a description of this WACCM run, which was conducted to 
support the Stratosphere-Troposphere Analyses of Regional Transport 2008 campaign. 

350 We specify the vertical distribution of Br y (the total bromine content of all inorganic bromine 
species) input to the PSS model, based on a relation between Br y and CFC-12. Profiles of CFC-12 
are from the NASA Global Modeling and Assimilation Office (GMAO) Goddard Earth Observing 
System Data Assimilation System Version 5 (GEOS-5) (Rienecker et al., 2007) assimilation for 
spring 2008. Salawitch et al. (2010) evaluated the sensitivity of stratospheric partial column BrO to 
355 the amount of stratospheric bromine supplied by VSL bromocarbons, termed VSL Br y . The value 
of VSL Br y is quite uncertain: Tables 1-9 of (WMO, 2011) gives a range of 1 to Sppt for this 
quantity. We use VSL Br y = 7 ppt because this value gives best overall agreement between OMI 
total column BrO and the sum of tropospheric partial column BrO inferred from ARCTAS and 
ARCPAC aircraft profiles and the simulated stratospheric BrO burden (see Auxiliary Material of 
360 Salawitch et al. (2010) for a discussion of the plausibility of VSL Br y = 7 ppt). 
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The PSS model is used to find BrO, at the local solar time of each OMI or GOME-2 orbit, as- 
suming production and loss of all species are in balance over a 24 h period of time for a stationary 
air parcel. Then, BrO abundance is integrated from the pressure of the tropopause (WMO definition 
of the thermal tropopause, based on GEOS-5 temperature) to 0.01 hPa, yielding stratospheric col- 
365 umn BrO. The diurnal variation of the stratospheric column BrO must be considered in estimating 
tropospheric BrO column, especially for the high SZAs of the level 2 data examined below. Fur- 
ther details of the PSS model are provided in Sect. 6.3.2 of SPARC CCMVal (2010), and references 
therein. 

As noted above, the Br y versus CFC- 12 relation used here consists of a baseline value of Br y plus 
370 7 ppt (representing stratospheric injection of bromine in the form of product gases of VSL bromo- 

carbons). Baseline Br y is set to zero at the tropopause and represents, above the tropopause, supply 
of bromine from the decomposition of methyl bromide (CH 3 Br), halons, as well as dibromethane 
(CH 2 Br 2 ). Other groups place CH 2 Br 2 into the definition of VSL Br y . We consider CH 2 Br 2 to 
be part of baseline Br y , however, because this compound is observed above the tropical tropopause 
375 (e.g. Wamsley et al., 1998). The total stratospheric Br y burden of 26 ppt used here is at the upper 

end of the present range of uncertainty (e.g. WMO, 2011). This formulation (7 ppt for VSL Br y plus 
baseline that includes CH 2 Br 2 ) was chosen because it results in best agreement between OMI total 
column BrO and the sum of stratospheric (modeled) and tropospheric (measured) partial column 
BrO. A wide range of values for VSL Br y yields “reasonable agreement” between these two terms 
380 when all uncertainties are considered (Fig. 6, Auxiliary Material, Salawitch et al., 2010). Despite 
the considerable uncertainty in stratospheric Br y and the impact of this parameter on tropospheric 
BrO column inferred from satellite observations of total column BrO, our approach of using a value 
for stratospheric Br y near the upper end is supported by the generally close quantitative agreement 
between inferences of tropospheric BrO column from the satellites and the in-situ data, described 
385 below. As noted previously, the fields of stratospheric BrO used here are approximately 27 % larger 
than those reported by Theys et al. (201 1) (with percentage differences that vary considerably with 
respect to location). 

2.4 Cloud parameters 

Optically thick clouds shield the underlying atmosphere from satellite sensors. The OMI rotational 
390 Raman (RR) cloud product (Vasilkov et al., 2008) is used to infer information about the shielding 
effects of clouds over snow and ice (Vasilkov et al., 2010). This product provides an estimate of 
the scene (combined cloud and surface) pressure over snow and ice surfaces. The Near-real-time 
SSM/I EASE-grid daily global Ice and snow concentration and Snow Extent (NISE) data set (Nolin 
et al., 1998) is used to identify snow and ice -covered pixels. When the difference between scene 
395 and terrain pressure (AP C ) exceeds 250 hPa we infer that clouds have led to significant shielding of 
tropospheric BrO from the satellite sensor (Vasilkov et al., 2010), as discussed below. 


13 



We use level 3 (gridded) cloud optical thickness and cloud-top pressure retrievals from the 
MO Derate-resolution Imaging Spectroradiometer (MODIS) (Platnick et al., 2003) on the NASA 
Aqua satellite as a secondary check on our detection of shielding clouds. Analysis of cloud effects 
400 on the satellite-derived tropospheric BrO columns is presented in Sect. 3.2. 

2.5 Meteorological data sets 

Tropopause heights are inferred along the flight tracks using tropopause pressure and geopotential 
height profiles from MERRA, a reanalysis based on GEOS-5 system (Rienecker et al., 2007) pro- 
vided by NASA Global Modeling and Assimilation Office (GMAO). These tropopause heights are 
405 used to integrate the aircraft BrO profiles and are given in Tables 1, 2a and 2b. MERRA also pro- 
vides meteorological parameters including sea level pressure, wind speed, and planetary boundary 
layer height. These meteorological data are used to infer relationships between BrO enhancements 
and the meteorological conditions. MERRA parameters are provided at 0.5° latitude x 0.667° longi- 
tude resolution. All parameters except geopotential height profiles are provided hourly; geopotential 
41 0 height is given every 6 h. 


3 Derivation of tropospheric BrO vertical column density 

3.1 Residual method to obtain tropospheric BrO vertical column density 


Tropospheric column BrO can be obtained from satellite total column BrO retrievals using the resid- 
ual method (e.g. Theys et al., 2011, and references therein). Here, as in Theys et al. (2011), we cal- 
41 5 culate tropospheric BrO vertical column densities (VCDx rop ) for each OM1 or GOME-2 pixel using 
the derived BrO total slant column (SCDxotaiX an estimate of stratospheric BrO vertical column 
(VCDstratX stratospheric air mass factors (AMFgtratX and an estimate of the tropospheric air mass 
factor (AMF'Trop) as follows: 


420 


425 


vrn SCDtiqp (SCD T otai-SCD Stfat ) (SCD T otai - VCD Str at • AMF strat ) 

Trop AMF Tf0 p AMFtvop AMFx r0 p ' 1 ' 

The SCD of a given absorber seen by a satellite sensor is defined as the amount of the absorber 
along an average light path taken by photons as they travel from the sun, through the atmosphere, 
and back to the sensor. The SCD is affected by scattering and absorption within the atmosphere as 
well as reflection off the surface and clouds. 

For a given altitude range (denoted by a subscript z ), The air mass factor (AMF 2 ) is used to 
convert SCD 2 to VCD 2 as follows: 


AMF 2 = SCD 2 /VCD 2 . (7) 

The sensitivity of UV radiance measurements to the BrO layer amounts varies with altitude. This 
variation depends on viewing geometry (mainly solar zenith angle, SZA, and view zenith angle. 
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VZA), surface albedo, cloud effects, and the vertical BrO profile. We must account for this varying 
430 sensitivity in the AMF. Using the optically thin absorber assumption, the AMF, can be formulated 
as 


AMF*=SCD*/VCD 2 = 


fW(z)N(z)dz 

jN(z)dz 


( 8 ) 


(Palmer et al., 2001 ; Theys et al., 201 1 ), where z is altitude, N(z) is the number density profile of the 
absorber, and W(z) is the weighting function profile that represents all the parameters influencing 
435 the AMF except the vertical profile of the absorber. 

Here, we use estimates of AMFstrat provided in the OMI and GOME-2 total BrO column prod- 
ucts. AMFstrat for GOME-2 is computed using the SCIATRAN radiative transfer model (Rozanov 
et al., 2005). The OMI algorithm uses a different approach to calculate AMFstrat; BrO VCDs are 
obtained by applying wavelength- and albedo-dependent AMFs to the BrO absorption cross sections 
440 prior to fitting. It is assumed that all BrO resides in the stratosphere in this procedure. An effective 
AMF is then defined as the ratio of the SCD (derived independently as described in Sect. 2.2) to this 
VCD. We assume that AMFstrat is represented by this effective air mass factor. 

We compute the tropospheric AMFxrop at 344.6 nm using the Linearized Discrete Ordinate Ra- 
diative Transfer (LIDORT) model (Spurr et al., 2001). We use a clear scene assumption in the AMF 
445 calculation. The DC -8 composite profile shown in Fig. 2a is taken as the default tropospheric BrO 
profile in the AMFxrop calculation. However, we use individual aircraft profiles of BrO, from both 
the DC -8 and WP-3D, to calculate AMFxrop lor point-to-point comparisons between aircraft in-situ 
and satellite-derived tropospheric columns. 

We have generated a look-up table of AMFxrop for various representative surface albedos and 
450 viewing geometries. A 5° interval for SZA and a 2.5° interval for VZA are used. Because the 
weighting functions are highly dependent on surface albedo, we use step- widths of 0.1 in the range 
[ 0 , 1 ] for surface albedo. 

The calculation of AMFxrop uses the derived OMI reflectivity at 33 1 nm for each pixel, from the 
OMI Total Ozone Mapping Spectrometer (TOMS) total O 3 product, as a proxy for surface albedo. 
455 We thus avoid using a surface albedo climatology, which can lead to AMF errors when the cli- 
matology differs from the actual surface albedo. Such errors can occur when the snow or sea ice 
distribution differs from the climatological mean, which is of particular concern at high latitudes 
during spring. 


3.2 Sensitivity of the derived tropospheric column to SZA, surface albedo, and clouds 

460 Here, we investigate the sensitivity of the satellite-derived tropospheric BrO retrieval to SZA, surface 
albedo, and clouds. Theys et al. (2011) analyzed the dependence of the weighting function on 
surface albedo and clouds for SZA — 45°. Here, we extend this analysis to a wider SZA range (40° - 
90°) relevant to polar observing conditions; at high latitudes (>60° N) during the day in the early 
spring SZA is >50°. Figure 3a shows that when the sun is relatively high in the sky (SZA <~60°), 
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465 OMI and GOME-2 should have good sensitivity to tropospheric BrO for surface albedos >0.5. 
Sensitivity is significantly reduced for darker surfaces (e.g. albedo — 0. 1). Sensitivity to tropospheric 
BrO decreases with increasing SZA; there is significantly lower sensitivity at 80° even for a surface 
albedo of 0.5. Figure 3b shows that at high surface albedo (0.9), there is increased sensitivity to BrO 
near the surface in addition to good overall tropospheric sensitivity for SZA up to 80°. 

470 The effect of optically thin clouds and aerosols on the sensitivity of OMI and GOME-2 to tro- 
pospheric BrO is dependent on surface reflectivity and viewing geometry. Vasilkov et al. (2010) 
showed that the sensitivity of UV satellite measurements to trace gas absorption near the sur- 
face in clear skies with moderately high surface albedo (70%) is approximately the same as 
for substantially cloudy conditions (optical thicknesses up to about 30) over a higher albedo 
475 surface (90%); the cloud shielding effect is much reduced for high albedo surfaces. In ad- 
dition, enhanced absorption takes place above a cloud of moderate to high optical thickness. 
According to Vasilkov et al. (2010), the UV cloud shielding effect is generally reduced over bright 
surfaces as compared to dark surfaces for nadir-viewing satellite sensors. For example, there is good 
sensitivity to tropospheric absorbers below thin clouds over bright surfaces near the nadir swath 
480 positions for SZAs <~67°. On the other hand, their simulation shows that there can be UV cloud 
shielding at the swath edge at high solar zenith angles, even for a thin cloud (r ~ 0.5). The clear 
scene assumption in the AMF calculation is appropriate over bright surfaces (i.e. snow/ice) when the 
derived scene pressure is close to the terrain pressure, which indicates that any particles along the 
line of sight must either be a thin, non-shielding cloud or an aerosol layer (the so-called Arctic haze) 
485 common in polar regions. The optical thickness of Arctic haze during ARCTAS was reported to be 
up to ~0.2 at 354 nm (Shinozuka et al., 201 1). Based on the calculations of Vasilkov et al. (2010), 
OMI and GOME-2 should have good sensitivity to surface BrO in Arctic haze conditions at the 
near-nadir swath positions, but may have reduced sensitivity at the swath edge. However, as shown 
by Theys et ai. (2011) and in Fig. 3b, optically thick clouds can shield satellite measurements 
490 from absorbers including BrO. 

Figure 4a shows computed tropospheric air mass factors. As shown for the weighting functions in 
Fig. 3, tropospheric air mass factors are higher over brighter surfaces owing to increased near-surface 
sensitivity. However, tropospheric sensitivity begins to be lost as SZA increases; tropospheric sen- 
sitivity almost disappears when SZA > 80° even for relatively bright surfaces. In this study, we 
495 present satellite-derived BrO tropospheric columns only for SZA < 80°. 

The sensitivity of AMFttop to the tropospheric BrO profile is assessed in Fig. 4b. For a bright 
surface (albedo = 0.9), no significant difference exists between tropospheric AMFs computed using 
the DC-8 composite profile and a profile where all tropospheric BrO is contained below 500 m 
for SZA <~60°. This is a consequence of good sensitivity at all tropospheric altitudes for bright 
500 surfaces and low SZA as shown in Fig. 3. For 60° < SZA <75°, the sensitivity to middle- and 
upper-tropospheric BrO is maintained while the near-surface sensitivity drops; this leads to a small 
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sensitivity to the profile shape under these conditions. Profile sensitivity decreases for higher SZAs 
as the sensitivity to the entire troposphere drops. For a darker surface (albedo = 0.4), the retrieval 
has lower sensitivity to BrO near the surface even when the sun is relatively high (SZA ~ 40°). 
505 In this case, the retrieval of total column BrO and our inference of the tropospheric column will 
be sensitive to the shape of the profile of BrO in the troposphere. If most of the tropospheric 
column happened to originate from BrO in the lowest 500 m of the atmosphere, then our 
inference of tropospheric column would be biased low over dark surfaces. 

We next examine tropospheric column BrO estimated from OMI in the context of the different 
510 sensitivities discussed above. Figure 5 shows OMI total, the model stratospheric, and the derived 
tropospheric BrO columns for OMI orbit 20050 on 22 April 2008 along with several parameters that 
impact the sensitivity of the radiance observation to the tropospheric BrO column. These parameters 
include 331 nm reflectivity (a proxy for surface albedo), AP C , and SZA. 

OMI total column BrO in Fig. 5a, within the region of the Barents Sea (black box), shows 
515 significant gradient. The region of high total BrO towards the upper left hand corner of 
the black box is associated with enhanced stratospheric burden (Fig. 5b, orange). Another 
region of large total BrO (lower portion of black box) is associated, by our analysis, with 
an enhanced tropospheric burden (Fig. 5c, red). This region of enhanced tropospheric BrO 
occurs over a bright portion of the Barents Sea (Fig. 5d, crimson, indicating snow or ice). 
520 Nearly zero tropospheric BrO column amounts are obtained over low surface albedo areas 
(OMI reflectivity < 0.5) of the Barents Sea (Fig. 5d, blue). Here, the retrieved total columns 
are generally less than retrieved columns over adjacent areas with higher reflectivity, leading 
to low tropospheric column BrO over parts of the Barent Sea. However, satellite-derived 
tropospheric BrO may not be reliable when the surface albedo is low. While it is possible 
525 tropospheric column BrO was truly low on 22 April over this region of the Barents Sea due 
to the lack of snow or ice leads that may be needed for bromine activation, it is also possible 
that our inference of low tropospheric BrO could either be a result of limited sensitivity over 
dark surfaces or an over-estimation of the stratospheric burden. The complicated sensitivity 
of satellite-derived tropospheric BrO to surface reflectivity and stratospheric burden requires 
530 concerted future study. 

In Fig. 5, the cloud shielding effect is shown over Siberia (marked with orange rectangles), where 
the OMI rotational Raman cloud product reports A P c > 250 hPa. Aqua MODIS also indicates cloud 
cover, with cloud top pressures between 450 and 500 hPa (not shown). Our derived tropospheric 
columns over this region are significantly less than those from adjacent areas without cloud cover. 
535 This supports the presence of tropospheric BrO mainly below 450-500 hPa and the ability of OMI 
and other nadir viewing sensors to capture information about tropospheric BrO over bright polar 
regions, with SZA up to ^80° when the sky is clear. Theys et al. (201 1) also used the cloud shielding 
effect to infer an amount of BrO in the troposphere, but applied this approach only in the tropics. 
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For completeness, we include all data in our analysis of aircraft measurements. However, we note 
540 aircraft soundings for which observing conditions likely produce large errors in the satellite-derived 
tropospheric BrO columns. In mapping surface BrO enhancement events with satellite retrievals, we 
use only data that has good sensitivity to the troposphere: SZA < 80°, A P c < 250 hPa, and surface 
albedo > 0.7. 

3.3 Uncertainty analysis 


545 For a single pixel, we express the uncertainty of the tropospheric VCD, ay cd-ivop’ * n a simplified 
form: 


^VCD T ,op 


( ^SCDTotaA 2 | { ^VCDstrat AMF g^at 
V AMFttop ) V AMFxrop 


(9) 


(Boersma et al., 2004; De Smedt et al., 2008; Theys et al., 2011), where crscD To tai * s uncerta i nt y 
of total slant column, <TscD S trat * s the uncertainty of stratospheric slant column, and cr amFttop the 

550 uncertainty of tropospheric air mass factor. This representation assumes that the different types of 
uncertainties are uncorrelated with each other. 

We take the fitting uncertainty derived from observed minus fitted radiances for each pixel as 
cr scD^otai • This assumes that the SCDx^tai error has a zero mean. The average fitting uncertainty 
at latitudes greater than 60° N is about 18%. Here, we do not consider systematic errors in 
555 SCDxotai t> r VCD Total* Systematic error will generally result in either a geographically uni- 
form over-estimate or under-estimate of total column BrO (VCDxotai)* There is synergy be- 
tween systematic error in VCDxotai and our prescription of the contribution of VSL species 
to stratospheric Br y . If subsequent analysis shows the estimates of VCDxotai BrO used here 
are biased high by a considerable margin, then clearly we must use a smaller contribution 
560 to stratospheric Br y to derive similar overall magnitude of tropospheric BrO. However, the 
geographic distribution of tropospheric BrO will not be strongly altered due to this synergy. 
An exploration of the systematic error in VCDxotai and the implication for tropospheric BrO 
will occur following analysis of a ground-based, OMI BrO validation campaign conducted in 
Fairbanks, Alaska during April 2011. 

565 We estimate the uncertainty of the stratospheric slant column (<TscD S trat) hy multiplying the un- 
certainty of the stratospheric column BrO by the stratospheric AMF. The uncertainty of the strato- 
spheric column BrO results from a root-sum-squares combination of three terms: the uncertainty 
in chemical kinetics that govern the BrO to Br y ratio, the uncertainty in the dynamics that govern 
CFC-12 (and hence Br y due to Source Gas Injection), and the uncertainty in VSL Br y . The uncer- 
570 tainty in chemical kinetics is evaluated by varying the rate constant of individual chemical reactions 
(including J values) by the 1 -sigma estimate of uncertainty given by Sander et al. (2006). The most 
important chemical term is BrO + NO 2 forming BrNO; 3 ; this rate constant is uncertain by about a 
factor of 2 at 220 K. The uncertainty due to dynamics is found by repeating calculations for plus 
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and minus 4 % variations in the abundance of CFC-12, because comparison to aircraft observations 
575 showed CFC-12 from GEOS-5 was accurate to within ±4% in the lower stratosphere (Salawitch 
et al., 2010). Finally, the uncertainty due to VSL Br y is set at ±27 %, which represents the mean 
difference in stratospheric column BrO resulting from our approach to handling this term compared 
to the approach of They s et al. (201 1). 

Here, we neglect the uncertainty in the tropospheric air mass factor as it is relatively small corn- 
580 pared with the other error sources. For example, Theys et al. (2011) describe sensitivity tests that 
show use of a single wavelength for the weighting function profiles leads to an error of less than 5 %. 
We find that the profile dependence produces error between 7 and 13 % for bright surfaces. 

4 Results and discussions 

4.1 Comparisons of aircraft in-situ measurements with satellite retrievals 

585 In this section, comparisons between aircraft and satellite inferences of tropospheric column BrO 
are shown for cases where agreement is good, agreement is poor but the cause for disagreement is 
understood, and agreement is poor and not well understood. These comparisons are performed to 
evaluate the veracity of the satellite-derived tropospheric column BrO product. We present a series 
of figures for various ARCTAS and ARCPAC flights. For each, we show maps of the satellite- 
590 derived total and tropospheric columns and time series of various aircraft measurements, collocated 
satellite-derived BrO data, and estimated stratospheric BrO amounts with vertical error bars derived 
using the formulation given in Sect. 3.3. Satellite-derived data are shown for all conditions (e.g. 
even in the presence of thick clouds that likely shield the surface). We also show plots of various 
parameters that affect the sensitivity of the satellite measurements to tropospheric BrO and note 
595 times and places where satellite observing conditions are not ideal for inferring reliable tropospheric 
BrO amounts. 

4.1.1 Cases of good agreement between satellite and in-situ data 

Figure 6a shows maps of the OMI and GOME-2 measurements of total column BrO as well as maps 
of tropospheric column BrO inferred from each satellite sensor on 17 April 2008. The track of the 
600 DC-8 and locations of BrO profiles, for the ARCTAS flight on this date, are marked. The bottom 
panel shows geographic regions of interest. OMI and GOME-2 orbits closest in time to the ARCTAS 
flight are used. 

The tropospheric column BrO maps have different spatial structure than maps of total column 
BrO. For example, total column BrO is highest over the north coast of Canada. These local maxima 
605 are significantly reduced in the tropospheric columns. A long tail of enhanced values near the North 
Pole is pronounced in the maps of tropospheric column BrO and is not as prominent in total column 
BrO. As discussed by Salawitch et al. (2010) and Theys et al. (201 1), inferring tropospheric BrO 
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from the total BrO column requires an accurate estimate of the stratospheric contribution to the 
total column, because variability induced by stratospheric dynamics is comparable in magnitude to 
61 0 variability induced by the surface release of bromine. Figure 6a is a perfect illustration of this point, 
which as noted above has been overlooked in many prior studies. Essentially all derivations of 
tropospheric column BrO prior to Salawitch et ah (2010) and Theys et ah (2011) relied on the 
use of zonally fixed stratospheric BrO, which, as first pointed out by Theys et ah (2009), will 
lead to large errors. 

61 5 Figure 6b shows time series plots of DC-8 flight data and collocated satellite measurements on 17 
April 2008. We indicate time in hours relative to the starting date of the flight; therefore times 
greater 24:00 refer to the following day. DC-8 data include aircraft altitude, in-situ O3, CO, BrO, 
active bromine and soluble bromide. Satellite data include OMI-derived BrO columns (total and 
tropospheric), A P c , OMI 331 nm reflectivity, and Aqua MODIS cloud optical thickness. Estimated 
620 errors of satellite-derived BrO columns are presented as vertical error bars. The tropopause height, 
as given in the MERRA data set, is at approximately ~ 7 km during the flight (not shown). CO and 
O3 data suggest that there were no stratospheric intrusions for any of the collected tropospheric BrO 
profiles. Obvious signatures of stratospheric air (low CO and high O3) were occasionally seen, but 
only when the aircraft was flying near or above the tropopause. 

625 Satellite-derived BrO tropospheric columns (orange line) are ~zero near Fairbanks (see top panel 
of Fig. 6b). As discussed in Sects. 3.1 and 3.2, tropospheric BrO information is not reliably retrieved 
when surface reflectivities are less than <~0.4, which is the case near Fairbanks (see the fourth panel 
of Fig. 6b). 

Both total and tropospheric column BrO show a slight dip near the pole, where the OMI and 
630 MODIS cloud products indicate shielding clouds (middle of Fig. 6b) . Here, the low values of derived 
tropospheric BrO over optically thick clouds likely result from the cloud shielding effect (Theys et 
al., 2011). Such clouds may also reduce tropospheric BrO columns by slowing photochemical 
production of active bromine. On the other hand, high values of satellite-derived tropospheric 
column are found along the DC-8 flight path both before the descent in altitude at ^24:15 
635 UTC and after the plane has ascended. Data from both OMI and MODIS indicate optically 
thick clouds were present during and after the ascent. The high value of satellite-derived 
tropospheric column during and after ascent could be due to the presence of BrO above or 
within these clouds. Profile #12 (see Fig. 2) shows a plume of enhanced BrO at about 3 km 
altitude, which could have been above the cloud deck. 

640 In-situ tropospheric column BrO, marked as black dots in the time series plot (the first panel of 
Fig. 6b), are obtained by integrating in-situ BrO mixing ratio profiles as described in Sect. 2.1. Un- 
certainties in the in-situ columns are shown as black vertical error bars. The horizontal bars indicate 
the distance covered by corresponding flight ascents/descents. OMI tropospheric BrO columns for 
point-to-point comparisons with aircraft data, shown as red squares, are obtained by averaging pixels 
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645 over the corresponding flight segment. 

The OMI tropospheric BrO columns exhibit a magnitude and variability similar to that of the 
DC-8 in-situ columns. For example, profile #12 (See Figs. 1, 2 and 6b) reports high BrO mixing 
ratios from the surface to km; the OMI tropospheric BrO column (as well as GOME-2, shown 
in Fig. 6 a but not in Fig. 6 b) is also relatively high at the same location. 

650 However, some discrepancies appear, perhaps because the aircraft captures small scale features 
that the satellite observations cannot resolve. For example, in-situ profile #13 (See Figs. 2 and 6b) 
shows a relatively large value for tropospheric column BrO, caused by enhanced BrO mixing ratios 
in a layer near 2 km. This could be a local event as the nearby aircraft column #14 does not show 
such an enhancement. The satellite retrieval agrees well with profile #14 but underestimates the 
655 abundance of BrO measured in profile #13. Overall, OMI tropospheric column BrO quantitatively 
agrees with the DC -8 in-situ column BrO to within the respective estimated uncertainties. The 
comparison between satellite and in-situ tropospheric column BrO for this flight is further quantified 
by Liao et al. (201 lb). These comparisons demonstrate the ability of OMI and GOME-2 to capture 
the magnitude and spatial distribution of tropospheric column BrO over bright surfaces for clear 
660 conditions. 

4.1.2 Examination of differences between satellite-derived and in-situ 
tropospheric columns of BrO 

We next describe two cases of apparent inconsistency between satellite and in-situ observations. 
Figures 7a and 7b show data corresponding to the DC -8 flight on 5 April 2008. In-situ measure- 
665 ments did not detect ozone depletion or high bromine over Hudson Bay, whereas OMI and GOME-2 
both report high total column BrO over vast regions of eastern Canada including Hudson Bay 
(Salawitch et al., 2010). BrO measurements were not available to produce profiles for this flight, but 
active bromine and soluble bromide were reported. Figure 7b shows that the derived tropospheric 
BrO column from OMI is low over Eastern Canada, where the total column is high. Similar results 
670 are found for GOME-2 (Fig. 7a). The OMI measurement of total column BrO shows very similar 
magnitude and structure as the calculated stratospheric column BrO. This suggests that the elevated 
total column BrO over Eastern Canada could be a consequence of high stratospheric columns. 
Furthermore, the in-situ CO and O 3 measurements indicate no evidence of stratospheric air below 
6 km altitude in ascending or descending flight tracks near Hudson Bay. 

675 Figures 8 a and 8 b show similar maps and time series, respectively, for the DC -8 flight on 8 April 
2008. GOME-2 observations during this flight are not available. Aircraft in-situ BrO measurements 
are also not available for this flight. This makes it difficult to evaluate the validity of the OMI-derived 
BrO tropospheric columns. The aircraft observed a severe ozone depletion event ([O 3 ] < 1 ppbv) 
and high levels of active bromine and soluble bromide near Alert (marked as black square), whereas 
680 the OMI total column BrO is relatively low and the inferred OMI tropospheric column BrO is at or 
near background levels. There are segments where inferred OMI tropospheric BrO is negative, 
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such as the early time period of 8. These values of negative BrO could be indicative of an over- 
correction for the stratospheric column, measurement noise for total coumn BrO, or some 
other error in calculation. Overall, we do not find extended regions of negative tropospheric 
685 BrO (either spatially or temporally). The ODE observed by the DC-8 on 8 April 2008 is thus not 
apparent in the satellite measurement of BrO on this date. 

As described in Eqs. (3), (4), and (5), BrO catalyzes the loss of ozone and, at the same time, 
BrO is produced by reaction of Br with 0 3 . This suggests that high BrO concentrations cannot be 
maintained if 0 3 is completely depleted, even though plenty of active bromine may be available. 
690 This explanation is consistent with other observations (e.g. Hausmann and Platt, 1994; Tuckermann 
et al., 1997; Neuman et ah, 2010; Friefi et al., 2011) that show high BrO concentrations only when 
ozone is partially depleted. The presence of a severe ozone depletion event is therefore indicative 
of an expectation of little or no enhancement of tropospheric BrO, despite the presence of large 
amounts of active bromine and soluble bromide. Liao et ah (201 la) also discuss the relationship 
695 between 0 3 and BrO; little BrO exists when 0 3 is low due to the partitioning between Br and BrO. 

OMI tropospheric column BrO shows enhancements in the vicinity of the 8 April 2008 ozone 
depletion event, during the 36 h time period prior to this flight (not shown). This enhancement is 
likely connected with the nearly complete removal of ozone observed by the DC-8 aircraft. 

4.1.3 Cases of disagreement between satellite and in-situ data 

700 Figure 9a shows maps of satellite-derived BrO columns and the flight track of the NOAA ARCPAC 
WP-3D on 19 April 2008. The OMI and GOME-2 data both indicate enhanced total and tropospheric 
columns of BrO over the Chukchi Sea. The aircraft flew into the area of elevated satellite-derived 
tropospheric column BrO at the northernmost extent of the track. 

Figure 9b shows time series plots for the flight of 19 April 2008. Two types of disagreement 
705 between in-situ and satellite tropospheric column BrO are observed. First, satellite-derived tropo- 
spheric column BrO is considerably less than the in-situ column over surfaces with moderate reflec- 
tivity in the presence of clouds. For example, the OMI estimate is less than the WP-3D columns (#19 
and #23) near Fairbanks. These observations were obtained in a region where the surface albedo is 
between 0.4 and 0,6, in the presence of low level clouds. This type of disagreement may be expected 
710 from the sensitivity studies described in Sect. 3.2. It is a challenge for satellite sensors to accurately 
retrieve tropospheric BrO over a moderately reflective surface in the presence of clouds. 

The second type of disagreement during this flight is illustrated by profile #21 (See Figs. 1, 2 
and 9b) , where both OMI and GOME-2 tropospheric column BrO are significantly enhanced (only 
OMI is shown in the first panel of Fig. 9b), whereas the in-situ BrO column does not show evidence 
715 of enhancement. Artifacts from clouds or surface albedo are not the likely explanation for this 
disagreement, as the surface was bright and no clouds are evident in A P c (fourth panel. Fig. 9b). 
The highest level of active bromine is observed from the aircraft sensor at the same location where 
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tropospheric column BrO from both satellite sensors show enhancements. O 3 mixing ratios of 
^lOppbv suggest the aircraft sampled a partial ozone depletion event (e.g, Ridley et al., 2003). 
720 Prior DOAS observations in the Arctic suggest that BrO should still be present for O 3 levels of 
~10 ppbv (Hausmann and Platt, 1994; Tuckermann et al., 1997). The disagreement between satellite 
and aircraft BrO could be related to the different spatial scales, and/or vertical coverage spanned by 
the respective instruments: the aircraft does not observe the complete profile and the field of view 
encompassed by the satellite covers a larger area than that sampled by the airplane. The aircraft 
725 may have missed an important part of the BrO profile as the minimum sampled height was 
151 m above the surface. Timing may also be a possible explanation for the disagreement, as the 
aircraft flew near the location of the enhancement about five hours after the OMI overpass and about 
three hours after that of GOME-2. 

Our analysis indicates that the large enhancement in total column BrO seen by OMI and 
730 GOME-2 over the Chukchi Sea on 19 April 2008 was tropospheric in origin, rather than strato- 
spheric, The WP-3D aircraft recorded highly elevated active bromine and partially depleted O 3 near 
the surface at this precise location, indicating recent association with elevated BrO. It is possible 
the satellite perturbation was caused by the presence of BrO at higher altitudes than those sampled 
by the aircraft at this location. By the time the aircraft reached 3 km altitude, it had left the region of 
735 highly elevated satellite BrO (terminal point of profile #23 (See Figs. 2 and 9b) ). Salawitch et al. 
(2010) presented an analysis of ground-based Max DOAS observations of BrO for April 2008 that, 
together with ARCTAS and ARCPAC BrO profiles, showed important contributions to column BrO 
often originate within the troposphere from altitudes above the top of the planetary boundary layer. 

Figure 10a and 10b shows maps of the BrO columns and time series plots for the NOAA WP- 
740 3D flight of 21 April 2008, respectively. The maps show enhanced OMI total column BrO over 
the north shore of Canada and Canadian Archipelago. The stratospheric model shows only a slight 
enhancement for these regions; the subtraction of the stratospheric model from the OMI column 
shows significant levels of tropospheric column BrO for much of the flight segment. Time series 
plots indicate that WP-3D in-situ columns #25, #26, and #27 (See Figs. 1, 2 and 10b) are lower than 
745 the OMI tropospheric column, but the differences are not significant given the uncertainties in OMI 
and in-situ tropospheric columns (i.e., error bars overlap). However, the absolute magnitude of 
the OMI tropospheric column is about a factor of 3 larger than the in-situ BrO column. There 
is, however, good agreement between in-situ column #28 and the satellite-based estimate. 

WP-3D profiles #26, #27, and #28 (See Figs. 1, 2 and 6b) were collected in close proximity (see 
750 Fig. 1), but show a large variation in the column amounts. As noted above, the profile #28 column 
agrees well with the satellite-based column whereas the other two in-situ based BrO columns are 
lower than the satellite-based estimates. As discussed above for other flights, this result may be 
explained by the fact that the aircraft captures small-scale spatial features while the satellite observes 
larger scales. 
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755 The overestimation of the tropospheric column BrO by OMI, compared to in-situ columns for 
profile #25, #26, and #27 (See Figs. 1 , 2 and 6b) , could potentially be explained by a stratospheric 
column that is not fully removed from the total column. One factor not considered in our analysis 
is irreversible, cross-tropopause exchange of of air parcels with elevated levels of Br y from the 
stratosphere to the troposphere (STE). Such transport events occur on the western flank of Arctic 
760 low pressure systems (Salawitch et al., 2010, and references therein). If the satellite signal were to 
originate from STE of Br y , the BrO signal associated with such air parcels may not have been 
sampled by the WP-3D because the maximum altitudes of the WP-3D profiles #25, #26 and #27 
(6 km) is well below the height of the tropopause (10.5 km). Future measurement of profiles of BrO 
from near the surface to the lowermost stratosphere, in the footprint of a satellite sensor after passage 
765 of an airmass with a low altitude tropopause, are needed to assess the importance of STE of bromine 
on the interpretation of the satellite record. 

4.2 BrO explosions observed with satellites 

Events of rapid enhancement of tropospheric column BrO are apparent from OMI and GOME-2 
observations after adjustment for the stratospheric burden of BrO, for the time period mid-March to 
770 late April 2008. The NASA DC -8 aircraft flew into an area of enhanced tropospheric column BrO 
near the North Pole on 17 April 2008 (Fig. 6 a) and the NOAA WP-3D aircraft flew near another 
tropospheric BrO enhancement on 19 April 2008 (Fig. 9a). Here we examine the event near the 
North Pole on 17 April 2008 and similar events in more detail using only OMI retrievals. The wide 
orbital swath and high spatial resolution of OMI, in addition to its frequent observations at high 
775 latitudes, provide a unique view of the temporal evolution of these events. 

Figure 1 1 shows the evolution of a tropospheric BrO enhancement event (“BrO explosion”) ob- 
served from 16 to 18 April 2008. Here, we only show observations when the following conditions 
are met, to provide reliable tropospheric BrO information as discussed in Sect. 3.2: SZA < 80°, 
reflectivity >0.7, and AP C < 250 hPa. The stratospheric column has been removed, as discussed in 
780 Sect. 3.1, using photochemical model output for the VSL Br y — 7 ppt simulation, for the local solar 
time of each OMI pixel. The major activation of BrO starts at ~22:00UTC on 16 April and lasts 
for ~30 h. We see activations near the North Pole and Canadian Archipelago. The spatial features 
of the elevated tropospheric column BrO change rapidly, with significant variations over the course 
of a day. 

785 We have ruled out a stratospheric origin for the enhanced total column BrO observations near 
the North Pole and Canadian Archipelago on these days. Since stratospheric BrO and O 3 columns 
exhibit a significant correlation (Theys et al., 2009; Salawitch et al., 2010), the OMI measurement 
of total O3 column can be used as a proxy for the spatial pattern of stratospheric column BrO. The 
comparison of the OMI tropospheric BrO column with the OMI total O3 column shown in Fig. 12 
790 indicates lack of correlation. Thus, the elevated region of total column BrO is not of stratospheric 
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origin. Figure 12 also shows few shielding clouds near the areas with high BrO columns. 

Figure 1 3 shows three tropospheric BrO enhancement events during April 2008 along with maps 
of sea level pressure, wind speed at 2 m altitude, and planetary boundary layer height from MERRA 
for the closest synoptic hour. For these events, the locations of enhanced BrO columns are co- 
795 incident with high near-surface wind speeds (3rd column) that are geostrophically consistent with 
localized low pressure systems. Neuman et al. (2010) also note high wind speed along the flight 
paths in conjunction with high concentrations of active bromine. Our analysis is consistent with 
the suggestion that strong surface winds associated with low pressure systems can trigger bromine 
activation via blowing snow (Yang et al., 2008, 2010; Jones et al., 2009; A. E. Jones et al., 2010). 
800 Figure 1 3 also shows that the spatial structure of high tropospheric column BrO is similar to that 
of the planetary boundary layer (PBL) height, although there is not always a precise alignment of 
these features. At high latitudes where the meteorological analysis is driven primarily by satellite 
data, the MERRA fields may contain displacement or other errors, particularly in near-surface fields. 
The following discussion should therefore be considered somewhat speculative in light of these 
805 uncertainties. We provide an animated visualization as a Supplement to depict evolution of enhanced 
tropospheric column BrO, sea level pressure, wind speed at 2 m and planetary boundary layer height 
from 16 to 18 April 2008. 

Neuman et al. (2010) reported a temperature inversion at 500 m during the BrO activation event 
of 19 April 2008, more or less consistent with the MERRA PBL height. Bromine chemistry may be 
81 0 related to PBL height as follows: BrO concentrations are controlled in large part by O 3 ; a sufficient 
amount of O 3 is necessary to maintain high BrO concentrations as explained in Sec. 4.1.2. Ozone 
is quickly consumed in a shallow boundary layer during bromine activation (Lehrer et al., 2004; 
Anderson and Neff, 2008). When O 3 is substantially depleted, production of BrO will cease and 
BrO will be destroyed by various reactions, including those in Eq. (4). On the other hand, the 
815 O 3 loss rate is quadratic in BrO, so a deeper boundary layer could slow down the chemical 
removal of 0 3 as BrO is diluted. 

The second row of Fig. 13 illustrates an event at 00:30 UTC on 18 April 2008 and the third row 
is for an event at 22:40 UTC on 19 April 2008. The two OMI orbits presented are closest in time 
to airplane flights into the elevated BrO and active bromine layers (profile #12 from DC -8 and #21 
820 from WP-3D, respectively (See Figs. 1, 2, 6 b and 9b) ). The locations of the profiles are shown as 
black diamonds on the maps of tropospheric BrO. The aircraft profiles of temperature, BrO mix- 
ing ratio, and active bromine mixing ratio are presented in the last column. Aircraft measurements 
show enhanced active bromine in the near-surface layer for both flights and enhanced amounts of 
BrO for the 18 April flight. Satellite-derived tropospheric BrO enhancements are closely related 
825 to near-surface parameters including sea level pressure, wind speed at 2 m, and PBL height. This 
observation suggests that the BrO activations originate at the surface, which of course is consistent 
with prior expectation. For the event on 18 April, a tail of enhanced BrO column exists parallel to 
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the wind direction inferred from sea level pressure. Begoin et al. (2010) suggested long range trans- 
port of tropospheric Br y for plumes in which recycling of bromine from condensed to gas phase 
830 sustains elevated BrO. The profile for BrO measured on 18 April by the DC -8 instrument (Fig. 13) 
also shows enhancements in the free troposphere, above the top of the PBL height (~500 m ). Si- 
multaneous profiles of CO and O 3 (not shown here, see Fig. 6 b) imply no stratospheric influence. 
The existence of elevated BrO above the top of the PBL could be due to vigorous convection over 
ice leads driven by warm exposed water, with BrO then dispersed horizontally by prevailing winds 
835 (Simpson et al., 2007b; Salawitch et al., 2010). Our results are consistent with a surface origin of 
elevated BrO as well as transport of Br y enriched air parcels away from the source of origin. In- 
deed, an animation of the panels in Fig. 1 1 (see Supplement) provides compelling evidence for this 
behavior. 

5 Conclusions 

840 We have estimated tropospheric column BrO for April 2008 from OMI and GOME-2 retrieved total 
column BrO retrievals and a model simulation of the stratospheric column BrO. The sensitivity of 
satellite radiances to tropospheric BrO depends on various parameters, such as surface albedo, solar 
zenith angle (SZA), viewing zenith angle (VZA), the shape of the BrO profile, and the presence of 
clouds. We have quantified these dependences with a radiative transfer model. A detailed compar- 
845 ison is provided between tropospheric column BrO found from OMI and GOME-2 measurements 
of total column BrO and tropospheric column BrO computed from in-situ observations of BrO and 
related species acquired by instruments aboard the DC -8 and WP-3D aircraft during the ARCTAS 
and ARCPAC field campaigns. 

Our analysis shows that retrievals of total column BrO from OMI and GOME-2, combined with a 
850 model estimate of stratospheric column BrO, can be used to retrieve realistic estimates of the mag- 
nitude and spatial variations in tropospheric column BrO provided that proper observing conditions 
are met. These conditions include bright surfaces (albedos >~0.7), low SZAs (<^80°), and the ab- 
sence of optically thick clouds. For the flight of 17 April 2008, satellite-derived tropospheric column 
BrO shows enhancements where the DC -8 reports intense bromine activation and moderate ozone 
855 depletion. The satellite and aircraft based measurements of tropospheric column BrO are in close 
quantitative agreement on 17 April 2008, particularly for DC -8 flight segments over bright surfaces 
in clear conditions. In contrast to prior expectation based on DO AS measurements, the tropospheric 
burden of BrO based on the aircraft profiles was distributed over a broad range of altitudes and was 
not restricted to the planetary boundary layer. The presence of elevated BrO above the PBL 
860 might be driven by convection over warm exposed water in regions of ice leads* with BrO then 
dispersed by prevailing winds (Simpson et al., 2007b; Salawitch et al., 2010). 

Many prior studies have used a zonally symmetric representation of stratospheric BrO to 
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derive tropspheric BrO from satellite observations. As put forth by Theys et al. ( 2009 ), Salaw- 
itch et al. ( 2010 ), and Theys et al. ( 2011 ) and reinforced here, this is a deeply flawed assumption. 

865 For the DC -8 flight of 5 April 2008, regions of enhanced total column BrO observed by OMI and 
GOME-2 near Hudson Bay are attributed to compression of stratospheric air associated with a low 
altitude tropopause, as neither the aircraft nor the satellite based measurements of tropospheric col- 
umn BrO show enhancements. The regions of enhanced tropospheric column BrO inferred from 
OMI and GOME-2 on 5 and 17 April 2008 are not readily apparent in maps of total column BrO 
870 because stratospheric variability imparts variations in column BrO comparable in magnitude to the 
tropospheric signal. Use of a zonally symmetric representation of stratospheric column BrO pre- 
vents the proper identification of the tropospheric BrO burden (e.g. Theys et al., 2009; Salawitch 
et al., 2010; Theys et al., 201 1). Caution should be applied when interpreting satellite-derived 
tropospheric BrO for SZA greater than ^80 degrees because, under these conditions satellite 
875 radiances have decreased sensitivity to absorption by tropospheric BrO (Fig. 4). 

Satellite and aircraft measurements of tropospheric column BrO do not always exhibit good agree- 
ment, at times for reasons that seem well understood and at other times for reasons that are unclear 
but may be related to differences in the timing or spatial coverage of the respective observations. For 
the severe ozone depletion event observed by the DC-8 near Alert on 8 April 2008, neither OMI to- 
880 tal column BrO column nor OMI tropospheric column BrO were elevated. Atmospheric conditions 
were favorable for the remote sensing from space of tropospheric BrO near Alert on this date (e.g. 
clear skies, high surface reflectivity). We believe the lack of a signal for OMI tropospheric column 
BrO near Alert on 8 April is consistent with our understanding of bromine chemistry: the production 
of BrO diminishes when ozone is severely depleted (e.g. Hausmann and Platt, 1994; Tuckermann et 
885 al., 1997). Hence, the association of elevated tropospheric column BrO and depleted O 3 is expected 
to be much stronger for a partial ozone depletion event (ODE) than a major ODE. Tropospheric 
column BrO from OMI does show an enhancement, near Alert, 36 h prior to the major ODE. The 
aircraft may have been capturing very small scale variability (several nearby profiles showed large 
differences) compared to space-based observations. Finally, aircraft observations generally did not 
890 sample the altitude region where stratosphere to troposphere transport of active bromine associ- 
ated with air parcels with elevated levels of inorganic bromine could potentially be affecting the 
satellite measurement of column BrO. 

We examined several events of rapid enhancement of tropospheric column BrO observed by OMI. 
Observations at high latitudes from this polar-orbiting, sun-synchronous satellite provide a unique il- 
895 lustration of the rapid temporal evolution of BrO activation because orbital periods are about 1 00 min 
apart. Our analysis indicates that BrO activation events occurring on 16, 17, and 18 April 2008 
are related to near-surface parameters such as low pressure systems, strong surface winds, and/or 
high planetary boundary layer. The satellite measurements reveal horizontal transport of activated 
bromine away from the source of origin, and the aircraft measurements show disbursement of BrO 
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900 within the free troposphere. The strong quantitative agreement between OMI and aircraft tropo- 
spheric column BrO on 17 April 2008 supports the validity of the rapid time evolution, on synoptic 
scales, revealed by the OMI tropospheric BrO product. The events on 17, 18, and 19 April suggest 
bromine activation via high winds over snow (Yang et al., 2008; Jones et al., 2009; Yang et al., 2010; 
A. E. Jones et al., 2010) as well as long range transport of Br y by surface winds (Begoin et al., 
905 2010). 

Supplementary material related to this article is available online at: 

http://\ @journalurI/\ @ pvoI/\ @fpage/\ @pyear/\@journalnameshortlower~\ @pvol-\ 

@fpage-\@pyear-supplement.zip. 
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Table 1. Information for NASA DC-8 BrO profiles shown in Fig. 1. “Day” means day of April 2008. Alts 
(AUt) means the bottom (top) altitude of the profile. “Col.” stands for BrO vertical column density (VCD). 
Units of BrO VCD are 10 13 molecules cm -2 . “Col. added'’ is the estimated amount of BrO that was added to 
the measured in-situ column to provide an estimate of the complete tropospheric column. 


No. 

Day 

UTC 

Lat. 

n 

Lon. 

t°] 

Alta 

[m] 

AltT 

[km] 

77trp 

[km] 

Col. 

IN-SITU 

Col. 

ADDED 

1 

16 

20:52 

69.7 

-145.3 

95 

4.2 

7.6 

4.36 

5.3% 

2 

16 

21:17 

71.1 

-147.4 

82 

4.2 

7.0 

2.85 

8.1% 

3 

16 

22:17 

73.2 

-159.0 

89 

4.0 

6.9 

3.97 

5.8% 

4 

16 

22:34 

72.1 

-161.9 

105 

3.6 

7.3 

2.14 

27.4% 

5 

16 

23:44 

66.3 

-165.5 

206 

4.1 

7.0 

3.04 

7.6% 

6 

16 

24:06 

65.2 

-165.2 

387 

4.2 

9.3 

2.87 

8.0% 

7 

17 

21:34 

70.3 

-148.7 

75 

5.2 

8.0 

2.68 

3.9% 

8 

17 

21:53 

71.0 

-148.2 

83 

6.3 

8.1 

2.18 

22.4% 

9 

17 

22:51 

77.1 

-148.4 

99 

6.3 

7.4 

3.86 

0.3% 

10 

17 

23:29 

80.3 

-148.4 

90 

5.5 

8.2 

3.44 

1.1% 

11 

17 

24:13 

84.9 

-148.5 

112 

5.4 

7.3 

3.04 

3.4% 

12 

17 

24:40 

87.2 

-148.4 

112 

5.9 

6.3 

4.35 

0.8% 

13 

17 

26:58 

78.3 

-156.5 

89 

7.0 

7.7 

5.20 

3.7% 

14 

17 

27:15 

77.1 

-156.6 

89 

4.5 

7.5 

2.57 

9.0% 

15 

17 

28:15 

71.8 

-157.5 

90 

4.1 

6.6 

2.11 

23.2% 

16 

17 

28:35 

71.3 

-156.8 

75 

5.8 

8.2 

1.88 

10.6% 
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Table 2a* Similar to Table 1 but for NOAA WP-3D profiles of BrO, measured on 12 to 15 April 2008. 


No. 

Day 

UTC 

Lat. 

n 

Lon. 

n 

AitB 

Cm] 

AltT 

[km] 

Htrp 

[km] 

Col. 

IN-SITU 

Col. 

ADDED 

1 

12 

21:14 

65.8 

— 148.1 

289 

4.6 

7.4 

2.55 

6.5% 

2 

12 

23:06 

70.6 

-152.7 

261 

5.2 

7.3 

1.31 

8.0% 

3 

12 

23:39 

72.6 

-153.3 

107 

4.1 

7.2 

1.96 

37.7 % 

4 

12 

24:00 

74.0 

-153.3 

77 

3.3 

7.0 

1.27 

43.8 % 

5 

12 

24:27 

75.1 

-151.2 

80 

5.0 

7.2 

3.34 

3.1% 

6 

12 

24:54 

75.1 

-144.5 

104 

5.0 

6.7 

0.15 

67.7% 

7 

12 

26:09 

74.0 

-139.7 

99 

3.4 

6.5 

1.41 

39.4% 

8 

15 

19:17 

64.8 

-149.5 

166 

3.5 

8.2 

0.62 

62.9% 

9 

15 

22:24 

66.4 

-165.0 

135 

4.6 

5.9 

0.73 

22.6% 

10 

15 

22:48 

67.8 

-165.0 

135 

4.3 

5.9 

0.81 

28.5% 

11 

15 

23:12 

69.1 

-165.1 

77 

3.2 

6.1 

1.71 

32.5 % 

12 

15 

24:30 

71.9 

-161.7 

93 

6.6 

6.7 

0.10 

2.9% 

13 

15 

25:00 

71.7 

-161.8 

100 

6.5 

6.6 

0.07 

17.9% 

14 

15 

25:45 

69.6 

-157.2 

207 

4.9 

6.3 

1.60 

15.2% 

15 

15 

26:50 

65.4 

-148.6 

203 

5.3 

7.3 

1.05 

10.0% 
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Table 2b. Similar to Table I but for NOAA WP-3D profiles of BrO, measured on 18 to 21 April 2008. 


No. 

Day 

UTC 

Lat. 

n 

Lon. 

n 

Alto 

[m] 

AHt 

[km] 

Htrp 

[km] 

Col. 

IN-SITU 

Col. 

ADDED 

16 

18 

22:39 

72.5 

-141.9 

92 

7.1 

8.5 

4.36 

-0.0% 

17 

18 

23:24 

73.0 

-137.7 

107 

6.8 

8.4 

0.85 

0.4% 

18 

18 

24:51 

70.9 

-144.2 

146 

5.2 

9.6 

1.92 

5.5% 

19 

19 

22:18 

65.5 

-148.5 

219 

4.8 

11.7 

4.49 

13.8% 

20 

19 

24:06 

71.5 

-156.5 

188 

4.7 

8.7 

2.44 

00 

vd 

21 

19 

26:51 

74.0 

- 160.0 

151 

3.8 

6.6 

1.79 

22.0% 

22 

19 

27:39 

72.6 

-153.8 

151 

7.4 

8.1 

2.12 

0.0% 

23 

19 

29:04 

66.7 

-149.1 

304 

5.4 

11.3 

2.29 

4.6% 

24 

21 

18:31 

66.0 

-145.6 

287 

6.1 

11.0 

4.38 

0.3% 

25 

21 

21:21 

72.7 

-127.3 

316 

6.3 

8.4 

1.43 

6.6% 

26 

21 

23:27 

69.6 

-136.7 

159 

4.6 

10.6 

0.96 

17.3% 

27 

21 

24:10 

69.5 

-136.7 

164 

4.7 

10.5 

1.31 

12.7% 

28 

21 

25:21 

69.3 

-137.6 

162 

5.9 

10.3 

2.97 

17.3% 

29 

21 

26:24 

66.1 

-146.0 

222 

5.4 

11.0 

1.58 

6.6% 



Table 3. Satellite-derived information corresponding to DC-8 BrO profile locations shown in Fig. 1. Units 
and abbreviations are as used in Tables 1, 2 and 2. “RatiooMi” means ratio of the OMI-derived tropospheric 
column to the in-situ column. “RatioGOME- 2 ” is the same but for GOME-2. 


No. 

Day 

SZA 

[°] 

R'iSlnm 

A P c 
[hPa] 

Col. 

OMI 

Col. 

GOME-2 

Ratio 

OMI 

Ratio 

GOME-2 

1 

16 

60 

0.80 

19.9 

0.98 

1.08 

0.23 

0.25 

2 

16 

61 

0.88 

46.7 

4.12 

3.91 

1.52 

1.44 

3 

16 

63 

0.91 

41.2 

2.99 

3.67 

0.78 

0.96 

4 

16 

61 

0,89 

44.2 

2.35 

2.78 

1.17 

1.39 

5 

16 

56 

0.85 

23.4 

2.15 

2.33 

0.74 

0.80 

6 

16 

55 

0.83 

49.3 

1.97 

2.26 

0.72 

0.82 

7 

17 

59 

0.85 

15.8 

2.70 

1.80 

1.03 

0.79 

8 

17 

60 

0.86 

21.3 

2.47 

2.65 

1.17 

1.26 

9 

17 

67 

0.88 

31.1 

3.10 

2.31 

0.81 

0.60 

10 

17 

70 

0.90 

46.3 

2.21 

1.54 

0.65 

0.45 

1 i 

17 

74 

0.87 

82.7 

3.11 

2.81 

1.05 

0.95 

12 

17 

77 

0.87 

146.9 

4.69 

5.08 

1.08 

1.17 

13 

17 

74 

0.90 

39.6 

2.85 

2.16 

0.55 

0.42 

14 

17 

75 

0.90 

32.9 

2.81 

2.05 

1.15 

0.84 

15 

17 

78 

0.86 

33.5 

2.48 

2.82 

1.26 

1.43 

16 

17 

80 

0.85 

38.6 

2.32 

2.17 

1.25 

1.16 
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Table 4a. Similar to Table 3 but for NOAA WP-3D profiles of BrO, measured on 12 to 15 April 2008. 


No. 

Day 

SZA 

[°] 

#331 urn 

A P c 
[hPa] 

Col. 

OMI 

Col. 

GOME-2 

Ratio 

OMI 

Ratio 

GOME-2 

1 

12 

57 

0.44 

11.9 

-0.46 

_ 

-0.18 

- 

2 

12 

62 

0.89 

85.9 

1.59 

- 

1.22 

- 

3 

12 

65 

0.89 

80.4 

1.24 

1.43 

0.63 

0.73 

4 

12 

67 

0.90 

89.5 

2.00 

2.31 

1.58 

1.82 

5 

12 

69 

0.90 

80.7 

1.58 

1.95 

0.47 

0.58 

6 

12 

72 

0.89 

70.7 

1.04 

2.00 

6.70 

12.94 

7 

12 

77 

0.90 

76.5 

1.23 

1.21 

0.87 

0.86 

8 

15 

62 

0.45 

155.9 

-1.39 

-0.42 

-2.23 

-0.67 

9 

15 

57 

0.85 

79.7 

2.09 

2.88 

2.85 

3.93 

10 

15 

58 

0.86 

92.9 

1.63 

2.20 

2.01 

2.71 

11 

15 

59 

0.86 

92.7 

2.19 

2.45 

1.28 

1.44 

12 

15 

64 

0.90 

104.3 

2.07 

3.56 

19.74 

34.07 

13 

15 

65 

0.91 

106.1 

2.33 

3.21 

32.63 

44.95 

14 

15 

67 

0.94 

106.7 

1.97 

2.22 

1.24 

1.39 

15 

15 

74 

0.46 

77.8 

0.28 

-0.58 

0.27 

-0.55 
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Table 4b. Similar to Table 3 but for NOAA WP-3D profiles of BrO, measured on 18 to 21 April 2008. 


No. 

Day 

SZA 

t°] 

Rs'Slmn 

A P c 
[hPa] 

Col. 

OMt 

Col. 

GOME-2 

Ratio 

OMI 

Ratio 

GOME-2 

16 

18 

62 

- 

- 

4.08 

2.71 

0.94 

0.62 

17 

18 

65 

- 

- 

4.20 

2.81 

4.97 

3.32 

18 

18 

67 

- 

- 

3.76 

2,14 

1.96 

1.12 

19 

19 

54 

0.61 

279.7 

0.91 

0.33 

0.20 

0.07 

20 

19 

62 

0.78 

176.6 

1.97 

2,55 

0.81 

1.04 

21 

19 

72 

0.83 

142.0 

4.23 

6.00 

2,37 

3.36 

22 

19 

77 

0.79 

179.7 

1.22 

1.80 

0.58 

0.85 

23 

19 

86 

0.72 

263.6 

0.78 

0.95 

0.34 

0.42 

24 

21 

63 

0.47 

266.0 

0.12 

0.25 

0.03 

0.06 

25 

21 

61 

0.83 

110.5 

3.03 

3.54 

2.13 

2,48 

26 

21 

62 

0.80 

198.1 

3.14 

3.29 

3.27 

3.43 

27 

21 

64 

0.80 

205.6 

3.02 

3.45 

2.32 

2.64 

28 

21 

69 

0.79 

193.5 

2.40 

2.87 

0.81 

0.97 

29 

21 

71 

0.47 

260.5 

0.12 

0.29 

0.07 

0.18 
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Fig. 1 . Flight tracks (lines) and the locations (numbers) of BrO profiles measured by instruments on board the 
(a) NASA DC-8 aircraft during ARCTAS and (b) NOAA WP-3D aircraft during ARCPAC. 


45 






Fig. 2. (a) Collection of DC-8 measurements of BrO (blue dots) and DC-8 composite BrO profile (black line) 
and (b) DC-8 BrO profiles #7, #10, and #1 1, (e) DC-8 BrO profiles #12, #13, and #14, (d) similar to (a) but 
for WP-3D measurements of BrO, (e) WP-3D BrO profiles #5, #22, and #23, (f) WP-3D BrO profiles #25, 
#26, #27, and #28. 
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Albedo = 0.1 & 0.5 Albedo = 0.9 




— A = 0.5, SZA = 40° 
A = 0.1, SZA = 40° 
A = 0.5, SZA = 60° 
A = 0.1, SZA = 60° 

A = 0.5, SZA = 80° 

— A = 0.1, SZA = 80° 



A = 0.9, SZA = 40° 

A = 0.9, SZA = 40° (cloud: 2.5km) 
A = 0.9, SZA = 60° 

A = 0.9, SZA = 80° 


Fig. 3. BrO weighting functions at 344.6 nm for various .solar zenith angles (SZA) and surface albedos (A), for 
(a) clear sky conditions and (b) clear and cloudy sky conditions, where the cloudy case is for an optically thick 
cloud (cloud extinction coefficient = 300) at 2.5 km altitude. 
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Solar Zenith Angle [deg] 


Solar Zenith Angle [deg] 


— Albedo * 1.0 — A=0.9, DC-8 composite profile 

— Albedo = 0.9 * A=0. 9, BrO mixing layer of 500m 

— Albedo = 0.8 — A=0.4, DC-8 composite profile 

Albedo « 0.7 A=0.4, BrO mixing layer of 500m 

•— Albedo = 0.6 
Albedo = 0.5 

— Albedo = 0.4 


Fig. 4. (a) Tropospheric air mass factors (AMFs) at 344.6 nm computed using the DC-8 composite BrO profile 
for various surface albedos as a function of SZA; (b) similar to (a) but also showing AMFs found assuming all 
of the BrO is below 500 m (dotted lines). A refers again to surface albedo. 
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Fig. 5. Measured total vertical column BrO, reflectivity, and cloud pressure for OMI orbit 20050 at 03:23 UTC 
on 22 April 2008 as well as calculated stratospheric column BrO for the time of OMI overpass and the inferred 
tropospheric residual column BrO. The panels highlight reduced tropospheric BrO over low albedo surfaces 
(black box) and in the presence of shielding clouds (orange box). See text for more details and discussion. 
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(a) April 17 2008 
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Fig. 6a. Maps on 17 April 2008 of (top left) OMI total vertical column density (VCD) of BrO, (top right) 
inferred OMI tropospheric VCD of BrO, (middle left) GOME-2 total VCD of BrO, and (middle right) inferred 
GOME-2 tropospheric VCD of BrO column. Negative values are included as dark blue. The flight track of 
the DC-8 on 17 April 2008 is indicated by black lines; aircraft profile locations where BrO was measured are 
marked with black diamonds. Geographic regions used in the paper are denoted in the bottom panel. 
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Fig. 6b. Time series plots for NASA ARCTAS DC-8 flight of 17 April 2008. For bottom three panels, line 
colors corresponding to the same color y axis. The top panel shows various estimates of VCD of BrO along the 
DC-8 flight track, including OMI total column, the modeled stratospheric column, and tropospheric columns 
based on aircraft in-situ measurements (black dots) as well as our analysis of OMI retrievals (red squares). The 
error bars of the tropospheric columns are described in the text. The two middle panels show various aircraft 
in-situ measurements along with the flight altitude. The bottom panel shows OMI reflectivity at 331 nm and 
cloud retrieved parameters from OMI and MODIS, as described in the text. Soluble bromide is not plotted 
when the reported mixing ratio was below the detection limit (1-6 pptv). 
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Fig. 7a. Similar to Fig. 6a but for the DC-8 flight on 5 April 2008. The GOME-2 measurement of total vertical 
column BrO is quite similar to OMI measurement and is therefore not shown. 
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Fig. 7b. Similar to Fig. 6b but for the DC-8 flight on 5 April 2008. Aircraft in-situ BrO measurements are not 
available for this flight because BrO was below the detection limit. 
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(a) April 08 2008 
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Fig. 8a. Similar to Fig. 6a but for the DC-8 flight on 8 April 2008. GOME-2 data are not available for this day. 
The location of severe ozone depletion and high active bromine sampled by the DC-8 is marked by the black 
squares. 
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Fig. 9b. Similar to Fig. 6b but for the WP-3D flight on 19 April 200S. Soluble bromide was not measured 
from the WP-3D. 







Fig. 10a. Similar to Fig. 7a but for the WP-3D flight of 21 April 2008. 
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Fig, 1 1 . OMI tropospheric BrO VCD for different orbits (UTC time as indicated) from 16 to 18 April 2008. 
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Fig, 12, (a) OMI tropospheric BrO VCD; (b) OMI total ozone column; (c) terrain minus cloud pressure for 
OMI orbit 19976 on 17 April 2008, 
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Fig. 13. First column; OMI tropospheric BrO VCD; Second column; MERRA sea level pressure; third col- 
umn: MERRA wind speed at 2 m altitude; fourth column: MERRA planetary boundary layer height; fifth 
column: aircraft in-situ profiles of BrO (red), active bromine (blue) and temperature (black). The rows corre- 
spond to OMI orbits. We show results for orbit 19979 on 17 April 2008 (top row), orbit 19990 on 18 April 2008 
(middle row), and orbit 20018 on 19 April 2008 (bottom row). Aircraft profiles were not obtained during 
the time that the satellite data shown in the top row were acquired. An animated GIF of this space-based 
component of this figure is available in Supplemental Material. 
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